Scattering theory of interface resistance in magnetic multilayers 
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The scattering theory of transport has to be applied with care in a diffuse environment. Here we 
discuss how the scattering matrices of heterointerfaces can be used to compute interface resistances 
of dirty magnetic multilayers. First principles calculations of these interface resistances agree well 
with experiments in the CPP (current perpendicular to the interface plane) configuration. 
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INTRODUCTION 



According to the scattering theory of transport developed by Landauer in 1957 JD , the conductance of a single-mode 



wave guide reads 



G -T— > W 



where T is the probability for a (spinless) electron at the chemical potential which approaches the sample from a 
reservoir on the left to be transmitted to a reservoir on the right. T = 1 represents perfect transmission and an infinite 
conductance. Much later, Buttiker realized that this expression should be replaced with 

G LB = (2) 

CJ ' 

When T — ► 1, the conductance is now finite, which physically represents the geometrical Sharvin resistance of the 
device. For the general multi-channel situation with spin the two-terminal Landauer-Buttiker formula reads: 

G LB = -j^ TjV.i* = \tja',ia\ 2 , (3) 

where tj a i ■ is the probability amplitude that an electron approaching the scattering region in transverse mode i and 
with spin a will be transmitted into an outgoing state with transverse mode j and spin a' on the other side of the 
sample. Some of the rather subtle points of this topic are discussed in [[|. 

The Landauer-Buttiker formula was introduced into the field of magnetoelectronics in order to compute the (mag- 
neto) conductance of microstructured magnetic multilayers ||. Asano et al. (| evaluated the effect of disorder 
numerically, whereas Brataas and Bauer |S used perturbation theory specifically for the Current Perpendicular 
to the interface Plane (CPP) geometry j| 0, ||, [)[ [l(|. The transmission coefficients are accessible to first-principles 
calculations, and G LB was computed for ballistic superlattices by Schep et al. p| , and for single, specular interfaces 
by Schep et al. p2|, van Hoof et al. [Of and Stiles and Penn pM. Recent improvements make it possible to extend 



^ , ab initio studies to disordered interfaces 15 



Schep et al. ]12|] pointed out that the point contact resistance R LB = 1/G LB is not the appropriate quantity to 
compare with the interface resistances, R 1 , which have been measured accurately by the Michigan State University 
■ collaboration for different systems 0. The experimental values are obtained by using the parameters of the series 
resistor model (in which interface resistances and bulk resistivities are simply added in series) to fit the measured data. 
The finding that the series resistor model works very well for the vast majority of samples is a strong indication for 
the absence of quantum size effects (see, however, the Discussion section), which can be rationalized by the disorder 
in metallic multilayers . In this "dirty" regime, transport can be described by semiclassical methods such as the 
Boltzmann equation which has been studied in detail for the CPP configuration by Valet and Fert |L6|. However, 
these authors did not take into account the discontinuity of the electronic structure at interfaces |ll[] which is the main 
source of the interface resistance jl2| . Schep et al. showed how the scattering matrix of an isolated interface should 
be incorporated into the Boltzmann equation and derived an explicit expression for the interface resistance in diffuse 
multilayers |1^] which can be directly compared with experimental results or numerical solutions of the Boltzmann 



equation 17 
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In the following we work out the details of the derivation of the interface resistance in a dirty environment, for 
which there was no room in the original article |T^ |, discuss the results in the light of new developments, and compare 
them with experiments. 

SEMICLASSICAL TRANSMISSION 

Let us begin by disregarding unnecessary complications associated with spin. At low temperatures, the distribution 
function for the CPP configuration in the presence of a weak electric field normal to a layered system and at a specified 
plane in a layer with index L can then be expanded as 



ft = fli + & fa* - e f) k - E F + t± 



(4) 



where i the state index, Ep is the Fermi energy of the total system and the superscript ± denotes whether the state 
i is right or left moving. /£ i is the equilibrium (Fermi-Dirac) distribution function. ji L — Ep is the local shift of the 
chemical potential, so /i L — fi L , is the potential drop between layers L and V . 7^ . describes the anisotropic part 
of the non-equilibrium distribution function, which vanishes in the sum over states. Once the 7^ , are known, the 
conductance can be calculated as 



G 



Total 



in 



(5) 



where Afj, Total is the potential drop over the total system. The spatially dependent distribution functions could be 
obtained from the general solutions of the linearized Boltzmann equation in the layers , determining the unknown 
linear coefficients by the boundary conditions at the interfaces. Here we will follow a different route by matching the 
distribution functions (0) at a given plane in each layer: 



(6) 



jeL> 



h.i - /J ( Rii ')y ft, j + ll ( T LL')ij /l', 
jeL jeL' 



(7) 



where the transmission (T) and reflection (R) probability matrices include bulk and interface scattering (see Fig. 1). 

In the absence of an applied bias all distribution functions are the same which leads to the current conservation 
condition: 



1 = £ (r LL ,) l3 + £ (R' LL ,) y - = £ (R M -) i3 + £ {T' LU ) t 

jeL jeL' jeL jeL' 



(8) 



and 
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jeL 



jeL' 



E( R "')«+E( T ^')i 

jeL jeL' 



Nr, 



N L , 



(9) 



(10) 



where Nl denotes the number of modes in layer L. Since the total transmittance does not depend on the sign of the 
applied bias: 



EE(t^- = EEw 

ieh' jeL ieL jeL' 

E ( R "') y ~ E ( R LL'h = N L -N L 



LL' 



(11) 

(12) 
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FIG. 1: 



Let us now locate L and L' at equivalent positions in different unit cells of a superlattice jT2|. The electronic 
structures are then identical and the scattering potential is symmetric, thus R — R' and T = T'. Substituting Eq. 
(i into (|0): 



ft 



\i = (ml - M E ( T ^'k + E (( T ™)ij ytj + &lv)h fi',j] 



(13) 



71,,: = (j*l> - Ml) E ( T ^)« + E (( R "' 

i j 

The superlattice symmetry implies: 



+ (Tll0 4J 7Z' j )- 



7L 



lL>,i = 7? 5 A M = Ml - Mi' 



which leads to 



(t). (J + (r; 
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We can now define a local conductance G s which can be written in a matrix notation as 



•7i 



A/Lt 



I-T + R) _1 T 



(14) 

(15) 
(16) 

(17) 



Note that the factor 2 is not related to spin. In the single mode limit this corresponds to Landauer's "old" formula 
(|j) , but differs from the many-channel generalization derived by Buttiker et al. |ll| . 

The transmission and reflection probability matrices of the superlattice unit cell in ( |l7j ) can be constructed via 
concatenation of the scattering probabilities of the individual constituents, i.e. interfaces and bulk materials. The 
transmission T12, T' 12 and reflection R12, R12 of a system of two scatterers (subscript 1 and 2) reads || [l{], [20| : 



Ti 2 = t 2 (i - r;r 2 ) 

R12 = Ri + TjR2 (I — R1R2) Ti 

T' 12 = T' 1 (I-R 2 R' 1 )" 1 T 2 

R' 12 = R 2 + T 2 R^ (I — R 2 R'jJ T 2 



(18) 
(19) 
(20) 
(21) 
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The semiclassical concatenation of matrices implies that scattered electrons lose phase memory when returning to 
the point of departure. As mentioned above, there is much evidence that magnetoelectronic devices to date are very 
dirty. Scattering at bulk layers and/or interfaces is therefore diffuse, which in three dimensional systems suppresses 
interference. Indeed it can be proven by Random Matrix Theory that transport is governed by the above concatenation 
rules |2l]]. This proof requires "isotropy" and is valid to first order in an expansion of 1/iV, where N is the number of 
conduction channels pjj . In transition metal systems this should be an excellent approximation, with the exception 
of few-atom point contacts or break junctions. 

The isotropy assumption of scattering means that the elements of the scattering matrix are equivalent and can be 
replaced by a constant value. If all elements of two matrices X (square, dimension N x N) and Y (dimension N x M) 
do not depend on their indices it is easy to show that 

Efo-XP'Y). = ■ (22) 

DIFFUSE TRANSPORT 

If the transmission though a bulk layer (material a) is ballistic we simply have (T Q )^ = Sij, (R a )y = 0. In the 
presence of diffuse scattering, the following choice parameterizes the material and layer thickness dependence: 

11 Is 

(T a ) tf = (T a ) y = j^YTZ'' {Ra) ^ = (Ra) « = J^YTZ (23) 

where 

e 2 pdN 



h A 



(24) 



consists of the bulk resistivity p, layer thickness d, and cross-section A. If we straightforwardly apply Eq. (||), we 
obtain for the dimcnsionless "point contact" conductance: 



lb _ h 



G L a B = E ( T «k = I \ Pada (25) 



ij N a h A 



where we recognize the (dimensionless) Sharvin resistance r„ = 1 /N a in series with the (dimensionless) conventional 
bulk resistance 

- (26) 



Substituting the expressions (g3|) into Eq. ( Jl T| ) , the conductance can be evaluated using Eq. (|22|) to yield the expected 
result: 

g s a = 2V[(I-T a +R a )- 1 Tj =— . (27) 



The transmission through an interface (I) -bulk (a) material is: 



T/ Q = T Q (I — RjRq,) x T/. (28) 

The matrix RjR a in the denominator stands for the multiple reflection of electrons between interface and bulk 
material. Ti a can be worked out under the random matrix assumption and 



1 , e 2 Pa d <* 

gf B h A 



2£f =£(T/«)y 

ij 

where now the Sharvin conductance N a in (|25| ) is replaced by the Landauer-Biittiker interface conductance: 



(29) 



9f B 

r 



= i = EI4f ( 3 °) 
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For a single interface a\I\(3 we need to evaluate 

T a ip — Tip (I — RoiRr^) 1 T Q (31) 
which proceeds along the same lines. The hetero point contact resistance is: 

rifp = ^-J^ v- =r a +r p + rf B . (32) 

The trilayer a|/|2/3|J|a has the transmission 

T a /2/3/a = T/jjq (I — R^j^R^/q) T a /^ (33) 

and resistance 

rinfiic = 2 h + r i B + r «) - r a h ~ 4 h (34) 
in contrast to the Landauer-Biittiker result: 

'' a/2/3/ 



LLi a = 2{r p + rf B +r a ) (35) 



In the series resistor model 0, 

AR T aI2 pi a = 2 (Afy + AR 7 + AR a ) (36) 
which agrees with the present equation if we identify: 



Ah ( r R r s a h + r s h \ Ah f 1 1 



i i 

+ a^ 



(37) 



which is the basic equation derived in |12|. The interface resistance in a diffuse environment is therefore not equal to 
the Landauer-Biittiker point contact resistance, but it can be obtained easily from it by substracting the geometrical 
Sharvin resistance! This result is readily extended to the two-channel resistor model, which holds for superlattices with 
collinear magnetization and sufficiently weak spin-flip scattering. Note that in this formulation the spin accumulation 
does not play a role. 

We can imagine a system in which bulk transmission is ballistic and interfaces are specular. Eq. ( |l7|) can then 
be calculated directly from the interface scattering matrix alone, two of which have to be concatenated in the case a 
single barrier [El| . That result would hold when the interfaces consist of specular islands much larger than the Fermi 
wave lengths, but thickness fluctuations of the layers which introduce the dephasing necessary for the semiclassical 
method H HI- 

DISCUSSION 

The method described here relies on several assumptions, the validity of which must be tested. From a theoretical 
point of view, it is satisfying that as far as the transmission probability matrix is concerned, our results agree with 
Random Matrix Theory Furthermore, it is important to note that the expression for the interface resistance 

( |37| ) does not depend on the bulk layer thickness or resistivity, in agreement with experiment. The assumptions 
of the model can be tested directly by comparison with numerically exact calculations for finite disordered systems 

are based on the surface Green's function method 



Because a minimal basis set is used, we are able 



J15[ . The parameter-free calculation of transport properties 15 
implemented with a tight-binding linear muffin tin orbital basis [ i 
to carry out calculations for large lateral supercells and model disorder very flexibly within these supercells without 
using any adjustable parameters. The electronic structure is determined self-consistently within the local spin density 
approximation. For disordered layers the potentials are calculated using the layer CPA approximation | p2[ . The 
interface roughness was chosen as a bilayer of 50 %/50 % alloy, but the results are not sensitive to moderate variations 
of alloy thickness and concentration, like a 60 %/40 % 40 %/60 % alloy concentration profile or between one and two 
mixed interface layers. 

Although the derivation of the interface resistance was based on the assumption of superlattice periodicity, this 
construction should only be viewed as a convenient model for a diffuse environment which is in fact less restrictive. 
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In a diffuse medium, the periodic boundary condition imposed above on the r y i can be replaced by the equality of 
the state-averaged JY 7i on both sides because the scattering matrix of a single interface sandwiched by diffuse bulk 
layers does not depend on the state index. It follows that in this limit Eq. ( |37| ) holds for a single interface as well. 

The isotropy assumption underlying the series resistor model appears to be hold for most structures. It is still 
interesting and important to explore their limitations. Recently, theoretical evidence has been obtained that the 
resistor model can break down [I?], |24j] . Deviations can occur in terms of a breakdown of the relaxation 

time approximation within the Boltzmann formalism E^j, or in terms of quantum corrections jl^, 15, ^4). These 



correction become significant for very clean or thin samples, but there is no consensus whether these have been 



observed in experiments on CPP spin valves |25, ^6[. Other complications may be residual spin-flip scattering at 
interfaces, which we cannot yet treat by first principles. 

Interface resistances are listed in Table 1 for a number of different systems. The differences between ARf and 
ARf B are very significant for highly transparent interfaces. The agreement of the computed interface resistances with 
experiments, which was already found to be good for specular interfaces of the Co/Cu systems is improved on 
including interface disorder jl5|. Recently CPP experiments on Fe/Cr |27| showed that the spin-averaged resistance 
agrees very well with the theoretical prediction, but that the polarization dependence is not yet completely understood. 

The model assumptions here differ from another semiclassical formalism, viz. the magnetoelectronic circuit theory 
p8|. There, the nodes in the circuit are taken to be at quasi-equilibrium (g^ i — 0), which is valid when the potentials 
drop exclusively over the resistive elements, which is the case for e.g. tunnel junctions, long wires, and point contacts, 
but not necessarily for magnetic multilayers. This implies, for example, that the angular magnetoresistance curves 
calculated by Huertas-Hernando et al. p9fl , should be modified along the lines expounded here 1 30 before a direct 
comparison with experiment [pl| can be carried out. 
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Table I: Results of first principles calculation of Eq. (|37|) jl5| and K. Xia et al., unpublished. 



system 


roughness 


ARi aj (mm 2 ) AR^jmin 2 ) 


AR^%(mm 2 ) 


AR£f„(fSlm 2 ) 


Au/Ag(lll) 


clean 




0.094 


0.094 


2.41 


2.41 


Au/Ag(lll) 


2 layers 50-50 alloy 


0.118 


0.118 


2.43 


2.43 


Au/Ag(lll) 


exp. 






0.100 ±0.008 


0.100 ±0.008 






Co/Cu(100) 


clean 




0.33 


1.79 


2.27 


3.11 


Co hcp /Cu(lll) 


clean 




0.60 


2.24 


2.67 


3.65 


Co/Cu(lll) 


clean 




0.39 


1.46 


2.39 


2.80 


Co/Cu(lll) 


2 layers 50-50 alloy 


0.41 


1.82±0.03 


2.40 


3.14 


Co/Cu(lll) 


exp. 
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0.26±0.06 


1.84±0.14 






Fe/Cr(100) 


clean 




2.82 


0.50 


3.51 


1.45 


Fe/Cr(100) 


2 layers 50-50 alloy 


0.99 


0.50 


1.68 


1.45 


Fe/Cr(110) 
Fe/Cr(110) 
Fe/Cr(110) 
Fe/Cr(110) 


clean 
clean Q 
2 layers 50-50 alloy 
exp. |p7| 


2.74 
2.11 
2.05 
2.7±0.4 


1.05 
0.81 
1.10 
0.5±0.2 


4.22 
3.53 


3.17 
3.22 



